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ABSTRACT 



We present observations with the new 1 1 GHz radiometer of the COSMOSOMAS ex- 
periment at the Teide Observatory (Tenerife). The sky region between 0° ^ RA ^ 360° 
and 26° ^ DEC ^ 49° (ca. 6500 square degrees) was observed with an angular resolution 
of 0.9°. Two orthogonal independent channels in the receiving system measured total power 
signals from linear polarizations with a 2 GHz bandwidth. Maps with an average sensitivity 
of 50 /iK per beam have been obtained for each channel. At high Galactic latitude (|6| > 30°) 
the 1 1 GHz data are found to contain the expected cosmic microwave background as well as 
extragalactic radiosources, galactic synchrotron and free-free emission, and a dust-correlated 
component which is very likely of galactic origin. At the angular scales allowed by the window 
function of the experiment, the dust-correlated component presents an amplitude AT ^ 9-13 
/xK while the CMB signal is of order 27 /iK. The spectral behaviour of the dust-correlated 
signal is examined in the light of previous COSMOSOMAS data at 13-17 GHz and WMAP 
data at 22-94 GHz in the same sky region. We detect a flattening in the spectral index of this 
signal below 20 GHz which rules out synchrotron radiation as being responsible for the emis- 
sion. This anomalous dust emission can be described by a combination of free-free emission 
and spinning dust models with a flux density peaking around 20 GHz. 

Key words: cosmology: observations - cosmic microwave background - Galactic anomalous 
emission - Galactic foregrounds 



1 INTRODUCTION 

Since 1984, the Teide Observatory has hosted experiments to ob- 
serve the cosmic microwave background at 10-33 GHz on angu- 
lar scales from degrees to arcminutes. Located at this observatory, 
the COSMOSOMAS experiment consists of two radiotelescopes 
designed to acquire sensitive observations of the microwave sky 
at frequencies in the range 10-17 GHz with an angular resolution 
of ~ 1 degree. It was conceived as a high resolution extension of 
the Tenerife experiment JDavies et al. 19961 IHancock et al. 19971 
[Gutierrez et al. 2000l l with the goal to provide 10-17 GHz ground- 
based data of angular resolution comparable to the lowest fre- 
quency channel of the Wilkinson Microwave Anisotropy Probe, 
WMAP ( Bennet et al. 2003al . Observations in this frequency range 
can contribute significantly to a better understanding of the fore- 
ground emission and to improve the limits imposed on cosmo- 



logical parameters by present and future CMB experiments. This 
appears particularly relevant for missions like Planclj^ which are 
planned to achieve a sensitivity of order 1 /^K per resolution ele- 
ment and infer cosmological parameters with accuracy better than 
1%. 

There is also an increasing interest in sky maps at few tens 
of GHz to elucidate the presence of a new microwave emission 
process different to the three classical Galactic foregrounds, 
namely, synchrotron, free-free and thermal dust emission. First 
detected as a dust-correlated signal in the COBE/DMR data 
(Kogut et al. 1996), this microwave emission was initially 
interpreted as due to free-free radiation. However, statistical 
detections of dust correlated signals by experiments like CBI 
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( IMasonetal. 200U . Saaskaton dde Oliveira-Costa et al. 19971 1, 
OVRO iLeitch et al. 1997t . the Tenerife experiment 
dde Oliveira-Costa et al. 19 99', Ide Oliveira-Costa et al. 2004t 
and also in the Green Bank Galactic Plane Survey (4.85 GHz) 
dPinkbeiner et al. 2004b strongly suggested that free-free emission 
cannot be the only underlying process. Draine and Lazarian 
(1998a,b) alternatively proposed electric dipole radiation from 
spinning particles and magnetic dipole emission as mechanisms 
behind this anomalous microwave emission. 

A dust correlated signal has also been detected in the WMAP 
data (Lagache 2003, Bennet et al. 2003a) but its interpretation is 
still a matter of debate. Bennet et al. 2003a argue that all fore- 
ground emission observed at WMAP frequencies can be explained 
in terms of the three classical foreground components. However, de 
Oliveira-Costa et al. (2004) show that the synchrotron template at 
22 GHz derived by the WMAP team under the previous assumption 
correlates with the 10 and 15 GHz Tenerife data in a manner in- 
compatible with the frequency behaviour of synchrotron emission. 
More recently, Fernandez-Cerezo et al. (2006) using high Galactic 
latitude observations at 13-17 GHz obtained with the first radiome- 
ter (C0SM015) of the COSMOSOMAS experiment (over an area 
of more than 3000 dcg^) find evidence that below 20 GHz the 
spectrum of the dust correlated signal turns over just as predicted 
by spinning dust emission. Unambiguous evidence for anomalous 
microwave emission in individual astronomical objects has been 
found by the COSMOSOMAS experiment (Watson et al. 2005) 
in the Perseus molecular complex and by CBl in the dark cloud 
LDN1622 (Casassus et al. 2006). 

In thispaper, we present new results of the COSMOSOMAS 
experimenrl obtained with the 11 GHz radiometer, hereafter re- 
ferred as COSMOll. These observations cover a sky region of 
6500 deg^ with an angular resolution of ~ 1° and sufficient sensi- 
tivity to explore further the behaviour of the anomalous microwave 
emission and other foregrounds at high Galactic latitudes. While we 
find cosmic microwave background and extragalagtic radiosources 
dominate the temperature fluctuations at 1 1 GHz in our beam scale, 
we also detect a highly significant dust-correlated signal at this fre- 
quency which is compared with previous results in the same region 
using COSMO 15 and WMAP data. 



Table 1. Summary of geometric configurations and observed declinations 
(complete in RA) in this work. Tilt has always been 5°. Units in degrees. 



Period 


a 


13 


DECmin 


DECmax 


01/11/03-12/05/04 


26.00 


64.50 


30.80 


50.80 


23/05/04-10/06/05 


29.50 


64.50 


23.80 


43.80 



good overall performance has resulted in similar sensitivities to 
COSMO 15 in a much shorter time scale. In this work, we present 
COSMOll and analyze data taken from November 2003 to 10**" 
June 2005. By stacking the best first 150 good observing days, 
COSMO 1 1 achieves a comparable sensitivity, 40-50 fiK per beam, 
to that of C0SM015 data (Fernandez-Cerezo et al. 2006b . 

The COSMO 1 1 radiometer was equipped with two similar re- 
ceivers, working at almost identical frequencies, but sensitive to 
orthogonalpolarizations in order to provide a better control of sys- 
tematics and explore potential polarization measurements. 

2.1 Experimental layout 

The optical and mechanical structure of the COSMO 1 1 experiment 
(see Fig.[TJ is very similar to that of C0SM015. The main differ- 
ences are the size of the circular, primary, flat, mirror of 3.0 m (in- 
stead of 2.5 m) and the aperture of the parabolic, secondary, mirror 
which is 2.4 m instead of 1.8 m. Both instruments were designed 
to achieve an angular resolution of 0.9°. A metal block with a 5° 
opening angle tilts the main, flat mirror. In this way as the mirror 
spins (at ~ 1 Hz), the observed path in the sky describes a circle of 
20° of diameter. During the period of observations relevant to this 
paper the inclination of the primary mirror has been changed twice 
while the position of the parabolic dish and the tilt of the primary 
mirror were kept fixed. A change of the inclination of the primary 
mirror causes a displacement in DEC of the center of the ring de- 
scribed by the beam on the sky. In Tab.[T]we give a summary of the 
values for the main geometrical parameters. The setup adopted for 
the observations led to data acquisition from declination 23.80° to 
50.80°. 



2 THE COSMOll RADIOMETER 

The COSMOSOMAS experiment consists of two independent sim- 
ilar radiotelescopes (COSMO 15 and COSMOll) aimed to obtain 
measurements of the microwave radiation in the frequency range 
10-17 GHz with average resolution of 0.9° (FWHM). The ex- 
periment is sited at Teide Observatory (2400m altitude, Tener- 
ife) and details of the observational strategy can be found in 
( [Gallegos et al. 200T] >. Essentially, each of the two instruments per- 
forms daily observations of a strip in the sky with diameter 20° in 
DEC, centred at the meridian. Earth rotation produces a map with 
complete RA coverage. An area equivalent to ~ 14% of the whole 
sky is then covered daily with a typical sensitivity of 1.0-2.0 mK 
per pixel, or roughly 0.3-0.7 mK per beam, after removal of most 
atmospheric contamination (see Sect. |3.1^ . 

COSMO 1 1 started scientific operation at the end of 2003, sev- 
eral years later than the higher frequency C0SM015, however its 

|http : / / www ■ iac_^es/pro ject / cmb/co5mosomas| of the Insti- 
tuto de Astroffsica de Canarias. 



2.2 Receivers 

The COSMOll receiver is a polarimeter, designed to receive or- 
thogonal linear polarisations through 2 total power channels (see 
Fig.|2j. 

The whole front-end of the receiver is housed in a cryostat 
and cooled to 15-20K by a dual-staged helium gas refrigerator. The 
phase centre of the feedhom coincides with the focal plane of the 
offset parabola. This provides the alignment of the cryostat. The 
cryostat window is a sandwich of thin Mylar film and high density 
polypropylene foam. Secondary calibration is provided by a tem- 
perature controlled broadband noise source which enters the beam 
through a side lobe. This coupling is obtained by placing the cal- 
ibration probe close to the cryostat window but not in the main 
beam. The feed horn is a conical corrugated design for maximum 
efficiency and minimum side lobes. The feed horn is followed by 
an orthomode junction. Both are designed by ERA technology Ltd 
and are matched to the Channel master 2.4m offset parabolic an- 
tenna. The two linear polar channels feed two similar 10-12GHz 
cryogenically cooled amplifiers from Berkshire Technologies Inc. 
The noise temperatures are close to lOK and the gain is 25dB min- 
imum. The complete total power channels form direct detection 
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chains whereby a matched square-law detector is used to detect 
the microwave signal without first mixing the signal down to an 
intermediate frequency. The receiver response just before entering 
the detectors is shown in Fig. [5] The mean, weighted, frequency 
defined as: 

/"max 
up{u)du, (1) 

where p{i') is the normalized (by setting it's integral to unity) linear 
spectral response of the filters, yields 10.89 GHz for channel 1 and 
10.87 GHz for channel 2. The overall 1/f knee frequency is about 
3Hz, which is faster than the spinning flat plate frequency of IHz. 
The 1/f is removed from the data stream by eliminating sufficient 
Fourier harmonics. See section 2.3. The rest of the receiving system 
is analogous to that of COSMO 15 ( IFemandez-Cerezo et al. 2006t . 

2.3 Receiver orientation and scanning strategy 

In this work, we only deal with data obtained for the primary con- 
figuration with channel 1 and 2 polarizations oriented North-South 
and East- West, respectively, see Tab. |2] In order to perform po- 
larization studies, it is necessary to rotate the receiver system and 
take data with a second orientation, e.g., with a 45° relative ro- 
tation with respect to the primary configuration. Observations in 
this second configuration are being currently conducted and po- 
larization measurements of the Galactic foregrounds and the cos- 
mic microwave background will be considered in a future paper. 



Table 2. Summary of the receiver configurations. The orientations referred 
to local meridian. Useful days are also given for Channel 1 and 2, respec- 
tively. 

Period Channel 1 Channel 2 Useful days 

01/11/03-10/06/05 North-South East- West 135,140 



although a faint polarized emission has been already measured in 
the Perseus molecular complex with the first data obtained so far 
dBattisteUi et al. 2 006). 

In order to minimize the effect of systematics in the data set, 
we also rotated the cryostat by 90°, which essentially produces an 
exchange in the receivers of channel 1 and 2. 

In Fig.[T8] the scanning strategy of COSMO 11 is shown. The 
circles in RA and DEC, the sky coverage due to Earth rotation, and 
the two polarizations, pointing towards the Celestial north pole: one 
North-South; the other. East- West. Since we are also interested in 
Galactic foreground emission, we also give a figure of the scanning 
strategy in Galactic coordinates. The sky image is that of the low- 
est frequency WMAP channel, i.e. K band. One can readily see the 
crossing of COSMOSOMAS rings with the Galactic plane, both 
center and anti center regions, and the high Galactic latitude region 
observed. We also note the position of Cygnus A, our main calibra- 
tor 
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We note that in 30 s the apparent motion of the sky corresponds to 
~ 0. 1 2° . Since the beam size is around 0.9° , there is no problem in 
averaging data every ~ 30 s. This is done mainly to save memory 
and allow faster data processing. Nevertheless, we store indepen- 
dently the data obtained for each spin of the mirror every ~ 1 s 
mainly for the technical monitoring of the data acquisition system. 
Every second, a diode signal of 2 K is introduced in the feed horn. 
This allows to measure changes in the receiver gain. Each day of 
data is saved into a FITS binary file. 

Every 30 s a stacked data set, called here a scan is stored. 
Each scan is comprised of 275 samples. In this way, each sample 
corresponds to an angular size of order nD/275, where D is the 
diameter of the circle subtended by the COSMO 1 1 beam on the sky. 
From Sect. |2.1[ D is ~ 20° and therefore the angular size of each 
sampling index is ~ 0.23°. Summarizing, the angular size of the 
data bin is of order 0. 12° x 0.23° , which is appropriate for the beam 
angular resolution of COSMO 1 1 (eventually, for the astronomical 
map we will choose a pixelization of 1/3° x 1/3°). 



3 DATA PROCESSING 

In this section, we describe the data analysis of COSMO 11 obser- 
vations, from the daily map making procedure and calibration to 
the cumulative maps corresponding to a period of approximately 
110 equivalent observing days for each channel. 



Figure 2. Receiver system of COSMO 1 1 . Two polarizations of total power 
signal at ~ 10.9 GHz are observed. 
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Figure 3. Relative spectral responses of the two tilter band passes in 
COSMO 11. 



2.4 Data acquisition 

The COSMO 11 data acquisition procedures are similar to those 
used in COSMO 15. Data are sampled over 4000 /is, leaving a 
blanking time of order 400 fj,s. A Fast Fourier Transformation is 
applied to the data for each mirror cycle, keeping 106 harmonics. 
In practice, we find that for the construction of the daily maps, it 
is enough to consider up to the 45**^ harmonic; higher harmonics 
do not contain astronomical information. Then, the data gathered 
every ^ 30s are stacked together, yielding a mean value for each 
of the 106 harmonics and their corresponding standard deviations. 



3.1 Daily processing 

The first step is the correction of gain fluctuations in the data. To 
that end, we divide the counts measured in each scan by the value of 
the 2 K calibration diode corresponding to that scan. A first clipping 
of bad data is performed on the scans consisting of rejection of 
those with saturated signals, no signal, spikes and/or glitches. The 
result is a sequence of scans with total power signal in acceptable 
ranges. 

In addition to instrumental noise, our data are affected by at- 
mospheric emission which causes a modulation in the total power 
every spin of the primary mirror (see Fig.[4|. After probing some 
techniques to remove this atmospheric modulation, we arrived to 
the conclusion that subtraction of a Fourier fit to each scan is a suit- 
able and efficient procedure. In this way, we also remove the main 
part of 1/ /-like noise. At 1 1 GHz, it suffices to remove some low 
order Fourier terms in order to get useful data (see later. Fig. |10[l. 
However, for the sake of a direct comparison with the previous 
COSMO 15 data, notably more affected by atmospheric emission 
than the current experiment, we choose here the same number of 
Fourier terms to be removed in each scan, i.e., a Fourier series fit 
of V"^ order (i.e., a constant term plus seven sin and cos terms). 
The fit of any scan is thus independent from the rest. In the Fourier 
series fit to each scan, we mask any single data sample within a ra- 
dius of 0.4 ° of a signal exceeding 4 times the standard deviation of 
the scan data. This is done to minimize the effect of bright sources 
on the fitted curve. After this individual scan cleaning a second 
clipping is applied consisting of rejection of scans with standard 
deviation 20 times larger than the mean standard deviation of the 
three preceding and three following scans. 

After this processing has been executed, we still find a much 
weaker daily modulation of the data taken at the same sky position 
(same line of sight or same sample index ). This modulation is pos- 
sibly associated to day-night atmospheric temperature variations. 
We correct it removing the first four Fourier terms in the fit to the 
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Figure 4. Typical scan obtained witli COSMO 1 1 . Up: atmospheric emis- 
sion appears cleariy as a long-scale modulation. The peak close to sampling 
index 85 corresponds to the observation of the Cygnus A region. Down: to- 
tal power after subtraction of the first seven Fourier harmonics. Structures 
to the left side of the Cygnus A region are associated to nearby galactic 
plane regions. 




Figure 5. Scan sequence for a typical day in the COSMO 1 1 data. Scans are 
cleaned. Notice the duplication of the observed sky region. The Galactic 
Cygnus region and Cygnus A are clearly observed near scan number 2000. 
The brightest sources seen along the first scans are the complex regions 
towards the Galactic anti-center 



whole day-night data set in any sky position. Here again, data in a 
0.4 ° interval centred in points which exceed the day-night mean 
value in 4a are masked in the Fourier series fit. The correction ap- 
plied to the day-night modulation modifies the window function of 
the experiment (? ). Because of the masking of the bright sources 
during the fitting procedure applied to the day-night modulation, 
the effect on point sources is very small, however CMB, very weak 
sources and diffuse foregrounds are affected and will be taken into 
account when generating adequate templates for cross-correlation 
with the experiment data. 



3.1.1 Daily map making 

Due to the scanning strategy of COSMOll, each astronomical 
source is observed twice (except in the upper and bottom sections 
of the scan). Fig.|5]shows the sequence of scans divided into its 275 
sampling sections (vertical axis) in terms of the scan number (time). 
The two point-like sources placed at approximately (2150,12) and 
(2250,85) correspond to the first and second crossing of Cygnus 
A with the beam of the experiment. Other sources are also seen 
daily, but with insufficient signal/noise ratio to be used as primary 
calibrators. For example, GB6 0319-1-4130 (3C84) can be seen at 
around (200,270). This source is known to show variability. From 
the 1980s until the epoch of the COSMOSOMAS observations, its 
flux has decreased from ~ 50 Jy to ~ 20 Jy. Our observations are 
in agreement with the flux extrapolated from the WMAP K band 
measurement (S = 11.1 ± O.lJy) and the spectral index published 
for this source by the WMAP team (a = -0.8 ± 0.1). The other 
bright sources seen in the map are complex structures in the Galac- 
tic plane towards the anti center. 

This double measurement of any sky region is then used in the 
construction of a daily map. First, knowing the local time of the 
entrance and exit of a known point-like source, say Cygnus A, we 
compute their hour angle (HA). Next, we use the values of the dec- 
lination of the source and also the direct measure of the inclination 
of the main, flat, mirror (a) and of the tilt angle (9) as starting val- 
ues to compute optimal values for the geometrical parameters that 
make to coincide the derived and observed HA. 
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Figure 6. Typical daily map. The Galactic Plane crossing is clearly seen, 
including the Cygnus A complex. Also clearly detected is the crossing to- 
wards the anti center region. 



Once, a, (3 and are inferred, we construct a daily map by 
projecting each of the observed pixels into a rectangular grid with 
RA in its horizontal axis and DEC in its vertical one (see Fig.|6]l. 
The chosen pixelization is 1/3° x 1/3°, which is adequate for the 
angular size of the sampling index and for the size of the beam. 
Due to the range of observed declination (20°), this method is ap- 
propriate and we do not resort here to more refined methods as the 
Healpix pixelization. 



3.1.2 Daily flux/temperature calibration 

Each daily set of data is calibrated in flux or antenna temperature 
with observations of the brightest stable point-like radio source in 
our maps: Cygnus A (see Fig.[6]l. In Fig.jTjone of the daily images 
of Cygnus A observed by COSMO 1 1 is plotted. 

We use the empirical model by dBaars et al. 19771 ) to derive 
the expected flux density (5*) from Cygnus A. For our frequency 
range it gives 



log Scyg [Jy] = a + b log //[MHz] , 



(2) 
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Cygnus A complex 



order 1.1-1.5 mK, close to the theoretical sensitivity limit of the 




Figure 7. Cygnus A as recorded in a single day of observation by 
COSMOll. 



where a = 7.161 ± 0.053 and b = -1.244 ± 0.014 which leads 
to a flux of Cygnus A of 5'cyg(10 - 12GHz) = 153 - 122Jy. 

For each daily observation of Cygnus A we check the beam 
of the instrument by fitting a two-dimensional Gaussian with an 
allowed base-label offset. In the frequency range of COSMOl 1, the 
Rayleigh-Jeans limit applies. Therefore, the anterma temperature 
associated with a given flux from Cygnus A is 



2k v'^Q. ' 



(3) 



where c is the vacuum speed of light, k is the Boltzmann's constant 
and n = (7r/41n2) FWHM^ is the solid angle subtended by the 
beam. 

The effect of the band-pass filters is to yield a weighted an- 
tenna temperature for each observation of Cygnus A according to 



r: 



CygA 



(4) 



is simply Eq. [s] 



where the script C stands for channel and Tcyg 
with S replaced by Eq.|2]and where p'~' [u) is the normalized linear 
spectral response of each filter Assuming a circular FWHM of 1 ° 
the antenna temperature results 1 10.2 mK and 1 10.9 mK for chan- 
nel 1 and 2 , respectively. For the fluxes, we get 137.9 Jy and 138.3 
Jy for channel 1 and 2, respectively. This shows that the differences 
due to non identical band-pass filters are of order ~ 0.5 % only. In 
Sect. |3.2| we plot the distribution of daily root mean square (RMS) 
in clean regions. 

On the other hand, for a beam size of 1° , the polarization of the 
Cygnus A region is less than 3-5 %. We have not attempted to per- 
form a detail study of the expected difference of the calibration for 
each of the two channels taking into account the COSMOSOMAS 
scanning strategy and data processing. The reason is that, as ex- 
plained later, the cumulative maps of each channel are recalibrated 
against unpolarized HII regions. Indeed, this results in less than 3 
% difference in the calibration of both channels for the Cygnus re- 
gion. As we will see, in this work we do not need to go to that 
precision in the data analysis, since, typical, statistical errors are of 
order 3-5% also. 

We remark that the error associated with the Baars et al. model 
is dominated by the rather large uncertainty in the spectral indexes 
a and b. However, we have checked that this model predicts flux 
values for Cygnus A at WMAP frequencies within 10 % of the 
WMAP observed fluxes. 

Finally, as an indication, for a good day, the mean RMS of 
data in regions placed at high Galactic latitude > 50°) is of 



3.2 Cumulative maps 

In this section, we present the first cumulative maps obtained with 
the COSMOll experiment. Maps are pixelized in a rectangular 
RA and Dec grid. Temperatures for each pixel are computed as 
the weighted contribution of all the observations relevant to that 
pixel obtained through the whole observing campaign. The mea- 
surements associated to each pixel are weighted each day with the 
inverse of the variance of the data obtained that day in a 3° x 3° 
area centred in the pixel under consideration. This applies to all pix- 
els except those possibly affected by relatively bright radio sources 
which are treated in a different way. In these cases, the daily weight 
is taken as the inverse of the variance of the data in high Galactic 
latitude regions where no strong radio sources are known to con- 
tribute. Pixels within the WMAP "KpO" mask (Bennet et al. 2003) 
are considered in this category. 

In Fig.[8]we present the final maps of COSMOl 1 obtained ac- 
cording to the data processing described in Sect. |3.1| i.e., removing 
a 7"^ order Fourier series along each scan and a 4*'' order one along 
day-night data to remove atmospheric contamination to the high- 
est possible level. As mentioned in the previous section, the map 
for channel 1 has the equivalent of 135 days of observations and 
for channel 2 of 140 days. The angular resolution slightly changes 
along this map. The values measured for Cygnus A in these two 
maps, a point source when compared with our beam, are listed in 
Tab. [3] In the final maps, for the region 120° < RA ^ 270° and 
25° ^ DEC ^45°, the RMS per beam (out of the masked re- 
gions) is approximately 52^K and 48jj,K for channel 1 and 2, re- 
spectively. 
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Figure 8. Stacked map of COSMO 1 1 channel 1 (top) and channel 2 (bottom) for the period November 2003 to June 2005. 
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Table 3. FWHM in stacked maps for Cygnus A 





FWHMra(°) 


FWHMdec(°) 


COSMO 11 channel 1 


0.80 ± 0.01 


1.03 ± 0.01 


COSMO 11 channel 2 


0.80 ± 0.01 


1.01 ±0.01 



Cygnus A complex (cumulative) 
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Figure 9. Top: Image of the Cygnus A complex for the period Novem- 
ber 2003 to June 2005. Bottom: Best 2-dimensional Gaussian fit to Cygnus 
A (dashed line). Observations (solid line) and difference with the fit (dot- 
dashed line). Cuts are plotted at the declination of the source. 



Finally, as mentioned in Sect. |3.1[ the COSMO 11 data are not 
so severely affected by the atmosphere as the C0SM015 data. This 
may allow a correction of the atmospheric containination using 
only low-order harmonics in the Fourier series and thus preserv- 
ing information on larger angular scales in the final map. In Fig.|10| 
we show the cumulative map obtained in a region of the Galactic 
Plane following the same fitting procedure described in Sect. |3.1| 
but restricted to use only first order harmonics in the Fourier series 
fit applied to each scan. The resulting maps (similar for both chan- 
nels) are clearly richer at large scales than the maps constructed 
with higher order harmonic subtraction. The full exploitation of the 
information provided by COSMO 11 at large scales is left for a fu- 
ture work. Here, we just verified that the values given in Tab. [3] 
are compatible with those measured in maps generated with atmo- 
spheric removal using different number of harmonics in the Fourier 
series. The result is rather independent on the number of harmonics 
which is taken as an indication that the implemented source mask- 
ing procedure works properly for our main calibrator. 
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Figure 10. Left: Galactic plane image in the cumulative map of channel 2 
where a order Fourier series has been subtracted along each scan. Right: 
Same, but with the present 7^^ order Fourier series atmospheric subtraction. 



3.3 Source recognition 

A search for radiosources was performed with SExtractor (Berlin 
& Amots 1996) on the map resulting from the combination 
of the two channels. Some of the brightest point sources are 
listed in Tab. |3.3| with its most likely identification. The coor- 
dinates given in the table are those provided by SExtractor and 
were used to identify the sources in the Kiihr (Ktih r et al. 198 U . 
3C (http://www.mrao.cam.ac.uk/) and Green Bank 
( [Gregory et al. 20()6| catalogues. We identify each radio source as 
the brightest source (> 1 Jy) found in these catalogues within 15 
arcminutes of the COSMOSOMAS pointing. The 1 1 GHz fluxes 
were obtained for these sources as in Fernandez-Cerezo et al (2006) 
by fitting the theoretical point spread function to the COSMO 11 
data. The errors listed in the table include the calibration uncer- 
tainty. The WMAP catalogue was also searched and 22 GHz fluxes 
are listed whenever available. The sources not found in the WMAP 
catalogue are either strong steep spectrum sources such as 3C123 
or weaker flat spectrum sources such as J2202-I-4216 which were 
too faint at the WMAP frequencies to be detected. 

Our results are generally consistent with independent 
observations of these radiosources at similar frequencies 
dKuhr et al. 198U , except for the source 3C84 for which we 
give much lower fluxes, but it is well known that the flux of this 
radiosource has progressively declined in the last few years. The 
flux measured by COSMO 11 is compatible with the tendency 
delineated by the observations at the UMRAO (Univ. of Michigan 
Radio Astronomy observatory) in the past decade. Many of the 
brightest radio sources show large variability and therefore the 
comparison of fluxes taken in different epochs can only provide a 
qualitative assessment of our data. 

We shall note that our ability to detect point radiosources in 
the COSMO 1 1 data at high Galactic latitudes is effectively limited 
to sources with fluxes above ~ 1 Jy due to statistical noise and 
to the "confusion noise" associated to an unresolved population of 
extragalactic radio sources. In Fig. [TT| a zoom of the final map of 
COSMO 11 shows the detection of several point-like sources near 
the Galactic plane in the anticenter direction. 

3.4 Noise analysis 

In this section, we describe the main features of the noise that is 
still present after subtraction of the first Fourier terms in the daily 
data. In Fig. [12] we plot the two-point correlation function for a 
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Table 4. A selection of the brightest sources detected in the COSMO 1 1 channels and the corresponding value in the WMAP K channel for comparison. 
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Source 


R.A. 


Dec. 


Channel 1 flux (Jy) 


Channel 2 flux(Jy) 


WMAPK 


3C84 


50.16 


41.54 


16.6 ± 1.7 


17.8 ± 1.8 


11.1 


3C123 


69.18 


29.33 


14.0 ± 1.4 


13.7 ± 1.4 




0923+39 


141.91 


38.89 


10.2 ± 1.0 


10.8 ± 1.1 


6.8 


J0418+3801 


64.64 


38.08 


9.5 ± 1.0 


9.0 ± 0.9 




3C345 


250.74 


39.86 


8.6 ± 0.9 


9.2 ± 1.0 


8.0 


J0555+3948 


89.17 


39.79 


3.9 ± 0.4 


4.7 ± 0.5 


1.2 


3C465 


354.73 


26.65 


3.2 ±0.4 


3.5 ± 0.4 
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Figure 11. A section of the final map of COSMOl 1 (channel 1) showing Figure 12. Two-point correlation function for a single day in COSMOl 1. 

the detection of point-like sources towards the Galactic anticenter. 



representative day. This is performed in a region RA=(1 15°, 267°), 
DEC=(25°, 45°) that corresponds to high Galactic latitude and us- 
ing the source mask described in Sect.|4] 

Some degree of correlation is expected between channels ob- 
serving simultaneously the same sky region. As mean covariance of 
the data of the two channels, we have obtained the average of the 
two-point correlation functions computed for each day whenever 
simultaneous observations in both channels were available (98 % 
of the days). The result is 0.40 mK, whereas the typical daily RMS 
of the data is 2.25 mK and 1.96 mK for channel 1 and 2, respec- 
tively. For comparison, the correlation between different days was 
found to be less than 1 % in good agreement with expectations for 
the amplitude of the astronomical signals. 



3.4.1 Stacked maps 

Finally, we have also considered the evolution of the noise per pixel 
when stacking daily maps. We verified that the RMS of the data 
contributing to each pixel decreased roughly as l/\/iV, where 
is the number of stacked daily maps. In the final maps of channel 
1 and 2, the typical errors in temperature (per pixel) are of order 
130-140 ^K. 



4 DATA ANALYSIS 

In this section, we study the degree of correlation of the previ- 
ous two COSMOl 1 stacked maps with multifrequency datasets 
available in the literature for the same sky region. The goal 
is to establish the contribution of CMB, synchrotron, brehm- 
stralung, extragalactic radiosources and anomalous microwave 
emission at the frequency and angular scales of our experiment. 
The set of maps to be correlated with COSMOl I are: (i) the 
3rd year maps for each of the five channels of the WMAP mis- 
sion i Hinshawetal. 2007( 5 

(ii) the extragalactic radiosource emis- 
sion maps constructed from the NVSS catalogue; (iii) the 408 
MHz map provided by (Haslam et al. 1982^ : (iv) the map pro- 
vided by (Finkbeiner 2003); and (v) the DIRBE maps at 100 ^m 
(DIRBE08) and 240 /im (DIRBEIO), and a HEALPix version of 
the dust map provided by i jSchlegel et al. I998| ) from LAMBDaQ 
at 100 /.im, (Aioo in the sequel). 

The method used for the computation of the correlation be- 
tween maps is described in detail in de Oliveira-Costa et al. (2004) 
and Fernandez-Cerezo et al. (2006)Correlation coefficients can be 
computed by comparison with just one template map or simultane- 
ously with several maps. The multi-correlation option is preferred 

^ These maps are given in the HEALPix pixelisation 
scheme iGorski et al. 20051 . 
^ http://lambda.gsfc.nasa.gov/. 
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Figure 13. Mask used for the cross-correlation studies, see text for details. 



when the templates present some degree of correlation themselves, 
as it is the case of those describing Galactic emission processes. 
In this work, all the correlation results except those corresponding 
to CMB come from a multi correlation analysis (unless otherwise 
stated) were the templates considered are: i) the V channel of the 
WMAP mission for CMB, ii) the NVSS map, iii) the 408 MHz 
map, iv) the Ha map and v) any of the three dust maps mentioned 
before. 

The correlation study considers data between 25.0° ^ 
DEC ^ 44.7° as allowed by the mask plotted in Fig. [TI] 
This is an extended version of the "KpO" mask used by WMAP 
where any adjacent pixel was also discarded for the correlation. 
Additionally, we masked some regions at high galactic latitude 
with strong radiosources to minimize any possible residual ef- 
fects of these sources after the adopted atmospheric filtering. The 
regions masked are (RA;DEC): (141.7° ± 3.7°, 38.7° ± 2.7°), 
(243.0° ± 3.7°, 34.3° ± 3.3°), (250.7° ± 3.7°, 33.5° ± 3.3°), 
(250.7° ±3.7°, 39.7° ±2.7°). 

Several other masking criteria were also considered in the cor- 
relation analysis: excluding all additional sources for which a di- 
rect extrapolation of the fluxes in the NVSS and GB6 catalogues 
yielded values higher than 1 Jy at 11 GHz; extending the KpO by 
two or three pixels around each masked source; and finally, exclud- 
ing any pixels with temperatures deviating several times the typical 
error per pixel of the map. The correlation results obtained with the 
various masks agreed within the statistical error bars of the method. 

In the following sections we deal with la statistical errors es- 
timated directly from the correlation method, but other sources of 
error should also be considered. Absolute calibration errors amount 
up to 10 % and chance alignment of regions with emission caused 
by physically unrelated processes also contribute to the final uncer- 
tainties. To determine the role of chance alignment, we have com- 
puted the correlation values and their dispersion for arbitrary rota- 
tions of the templates around an axis perpendicular to the Galactic 
plane. The dispersions of the correlation values were found usually 
1-2 times higher than the statistical error given by the correlation 
method. Only those correlation values larger than 2 times the statis- 
cal error will be considered a reliable result. 

The COSMO 11 maps are convolved to a final angular res- 
olution of 1.12° X 1.12° to set on equal foot the COSMOU and 
C0SM015 data, so both experiments can be compared directly. 
We have not considered the different polarization of the two 
COSMO 11 channels in producing the set of templates. This may 
lead to systematic differences in the correlation results obtained for 
each channel with the same templates. In general, these differences 



will not be larger than a few percent of the amplitude of the corre- 
lated signals, i.e. smaller than the statistical error of the correlation 
method. 

With the adopted atmospheric filtering, the window function 
describing the data of the COSMO 11 radiometer is analogous to 
that published by Femndez-Cerezo et al. (2006) for COSM015. 
We remark that the filtering applied to the COSMOSOMAS data 
during the reduction process is also applied to any other map used 
in the correlation analysis. This erases information on the largest 
scales, so the results need to be understood as the correlation level 
at the remaining scales, which essentially range from ~ 1° to 5°. 
Finally, here we have reprocessed the COSMO 15 data obtained 
by Fernandez-Cerezo et al. (2006), according to the procedures 
explained in Sect |2.4| which are slightly different to those used 
in the former work. The cross-correlation results obtained for the 
COSMO 15 data are also slightly different. 



4.1 Correlations with WMAP data 

InTab.H] we list the correlation values obtained for the COSMOSO- 
MAS (labelled as CI li and CI b for the two COSMOU channels, 
and C13, C15 and C16 for the three COSMO 15 channels) and the 
WMAP channels for three different galactic latitude cuts. First, we 
note that the correlations among the highest frequency channels of 
WMAP indicate the level of the CMB signal left by the window 
function associated to the adopted COSMOSOMAS data process- 
ing scheme. The correlations between the Q, V and W channels 
show a detection of CMB with amplitude 26-26.5 /iK at the vari- 
ous Galactic latitude cuts. This is in agreement with the expected 
CMB fluctuation after removal of the Fourier terms described in 
Sect. 13.11 The autocorrelation values for these channels are ~0.5 
higher as expected for the noise contribution. 

Approximately the same level of correlation (26-29 fiK) is 
found between the five COSMOSOMAS channels, and the high 
frequency WMAP channels (see Fig. |14^ . The CMB signal is de- 
tected in both COSMOU channels with a high significance (> 
25cr). CMB polarization is expected to cause differences between 
the correlation values obtained for each of the two COSMOU 
channels smaller than the statistical uncertainties. Therefore no at- 
tempt was made to measure this polarization effect in our data. 
When we consider regions at low Galactic latitude, as e.g. j6j > 
30°, the amplitude of the CMB correlated signal slightly decreases, 
this is an artifact due to the strong influence of the Galactic plane 
signal which affects the processing of COSMOSOMAS data. 

The amplitude of the correlated signal between COSMOU 
and the low-frequency WMAP channels increases to ~43 and 35 
^K, for the K and Ka channels (|6| > 40°), respectively. This is 
attributed to signals in common other than CMB. At high galactic 
latitude and considering the angular scales relevant to the present 
analysis, the expected major foreground contamination in com- 
mon with these WMAP channels are unresolved extragalactic ra- 
diosources as we will see in the next sections. 



4.2 Correlations with unresolved radio sources 

In order to determine the level of the contribution of unresolved 
radio sources in the COSMOSOMAS data we proceed as in 
Fernandez-Cerezo et al. (2006).Essentially, for each galactic lati- 
tude and at each frequency the correlation values of COSMOSO- 
MAS with the WMAP K map (values of Tab.|5j are corrected from 
the CMB contribution (obtained as mentioned above), as well as 



COSMOSOMAS 10-12 GHz Data. 1 1 



Table 5. Values of the temperature correlation in /^K for the COSMO 1 1 cumulative map with respect to the maps of the five WMAP channels for different 
cuts in Galactic latitude, b. Errors are 1 a. 



Template 


Clli 


Cll2 


CI3 




C15 




C16 


WMAP_K 




WMAP_Ka 




WMAP_Q 




WMAP_V 




WMAP_W 






|b|> 30° 


WMAP.K 


42-5 ± 0.9 


43 .8 ± .8 


32.1 ± 1.0 


30 


9 ± 1 


2 


31.9 ± 3.0 


29-3 ± 


1 


26.0 ± 


1 


25 


3 ± 


1 


24.7 ± 


1 


24.4 ± 


1 


WMAP_Kii 


33-0 it 0.9 


34.0 ± 0.8 


27.3 ± 1.0 


26 


3 ± 1 


2 


28.4 ± 3.0 


28-0 ± 


1 


27.4 ± 


1 


26 


5 ± 


1 


26-2 ± 


1 


26.0 ± 


1 


WMAP_Q 


29-9 ± 0.9 


31.1 ± 0.8 


26.1 ± 1.0 


25 


5 ± 1 


2 


27.7 ± 3.0 


27-5 ± 


1 


26 


8 ± 


1 


27.0 ± 


1 


26.4 ± 


1 


26.2 ± 


1 


WMAP_V 


27-0 ± 0.9 


27.7 ± 0.8 


23.9 ± 1.0 


23 


9 ± 1 


2 


25 .0 ± 3.0 


26.9 ± 


1 


26 


5 ± 


1 


26 


4 ± 


1 


26.9 ± 


1 


26.2 ± 


1 


WMAP_W 


25-1 ± 0.9 


25.7 ± 0.8 


23.0 ± 1.0 


23 


5 ± 1 


2 


24.9 ± 3.0 


26.5 ± 0.1 


26 


3 ± 


1 


26 


2 ± 


1 


26.2 ± 


1 


27.0 ± 0.1 


|fa|> 40° 


WMAP_K 


42-7 ±1.0 


44 .1 ± .9 


34.2 ± 1.1 


32 


8 ± 1 


4 


33 .3 ± 3.5 


29-1 ± 


2 


25 


8 ± 


2 


25 


2 ± 


1 


24-7 ± 


1 


24.4 ± 


1 


WMAP.Ka 


34-4 ±1.0 


35.6 ± 0.9 


29.9 ± 1.1 


28 


3 ± 1 


4 


29.9 ± 3.5 


27-7 ± 


2 


27 


1 ± 


2 


26 


2 ± 


1 


26-0 ± 


1 


25 .8 ± 


1 


WMAP.Q 


31-7 ±1.0 


33.1 ± 0.9 


28.8 ± 1.1 


28 


± 1 


4 


29.6 ± 3.5 


27-4 ± 


2 


26 


5 ± 


2 


26 


7 ± 


1 


26-3 ± 


1 


26.1 ± 


1 


WMAP.V 


29-0 ± 1.0 


29.7 ±0.9 


26.9 ± 1.1 


26 


5 ± 1 


4 


27.2 ± 3.5 


26-8 ± 


2 


26 


3 ± 


2 


26 


2 ± 


1 


26-8 ± 


1 


26.1 ± 


1 


WMAP.W 


27-3 ± 1.0 


28.1 ± 0.9 


26.1 ± 1.1 


26 


± 1 


4 


27.3 ± 3.5 


26-4 ± 


2 


26 


1 ± 


2 


26 


± 


1 


26-1 ± 


1 


26.8 ± 


1 


|fa|> 50° 


WMAP.K 


41-1 ± 1.1 


41.3 ± 1.0 


33.6 ± 1.3 


31 


± 1 


6 


33 .5 ± 3.9 


29-3 ± 


2 


26 


3 ± 


2 


25 


8 ± 


1 


25.3 ± 


1 


25 .0 ± 


1 


WMAP.Ka 


34-2 ± 1.1 


34.3 ± 1.0 


29.9 ± 1.3 


26 


9 ± 1 


6 


30.6 ± 3.9 


28-1 ± 


2 


27.4 ± 


2 


26 


7 ± 


1 


26.5 ± 


1 


26.2 ± 


1 


WMAP.Q 


31-6 ± 1.1 


32.1 ± 1.0 


28.9 ± 1.3 


26 


7 ± 1 


6 


30.7 ± 3.9 


27-8 ± 


2 


26 


9 ± 


2 


27 


1 ± 


1 


26.6 ± 


1 


26.4 ± 


1 


WMAP_V 


29-1 ± 1.1 


29.0 ± 1.0 


27.2 ± 1.3 


25 


3 ± 1 


6 


28.4 ± 3.9 


27-3 ± 


2 


26 


7 ± 


2 


26 


7 ± 


1 


27-1 ± 


1 


26.4 ± 


1 


WMAP_W 


28-0 ± 1.1 


27.9 ± 1.0 


26.9 ± 1.3 


25 


7 ± 1 


6 


29.0 ± 3.9 


26-9 ± 


2 


26 


5 ± 


2 


26 


5 ± 


1 


26-5 ± 


1 


27.1 ± 


1 
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Figure 14. CoiTelation of COSMOSOMAS and WMAP channels with the 
V channel of WMAP. EiTor bars are 1 a statistical errors. 



Figure 15. Predicted (circles) and observed (diamonds) contribution from 
unresolved radio sources, after taking into account COSMOSOMAS ob- 
serving strategy and data processing. See text, and Tab.|6]for details. 



from the corresponding correlation with the Galactic synchrotron 
(408 MHz map), free-free ( Ha template) and dust (DIRBE08) (see 
next sections for details on these correlations and their statistical 
errors). The correction consists in a subtraction of each of the pre- 
vious temperature correlation values in quadrature. We obtain the 
results listed in Tab. |6] and plotted in Fig. [TS] The data point at 23 
GHz in Fig.[T5]is obtained correcting the autocorrelation value of 
the WMAP K for the contribution of the CMB and the other fore- 
grounds. It is clear the existence of an additional signal in common 
between the COSMOSOMAS data and WMAP K band whose am- 
plitude is independent of galactic latitude and decreases with in- 
creasing frequency. At the frequencies of COSMOSOMAS and at 
the angular scales under consideration, this additional signal turns 
out to be the second most relevant contribution (after CMB) to the 
temperature fluctuations at high Galactic latitudes. 



In order to interpret these results we use the models by de Zotti 
et al. (2005)and estimate the expected contribution of unresolved 
radio sources. Assuming a detection limit of 1 Jy for resolved ra- 
diosources (masked in our correlation analysis) and considering our 
window function these models predict at 11 GHz a total contribu- 
tion of unresolved radiosources of order 29.0 fiK (Gonzalez-Nuevo 
private communication) For the COSMO 15 frequencies the mod- 
els predict contributions of order 21.0 /iK, 18.2 /iK, 16.3 /iK for 
C13, C15, and C17 channels, respectively. Finally, for the WMAP 
K channel, the predicted contribution is of order 1 1.3 /iK. The ex- 
cess temperature correlation values in Tab.|6]agree remarkably well 
with these model predictions. Both are compared in Fig.|15| 
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Table 6. Temperature cross congelation excess between the COSMOSOMAS and WMAP K map after contributions from CMB and Galactic emission are 
accounted, see text for details. All the maps are processed according to COSMOSOMAS strategy. Errors are 1 a. The error bars are computed taking into 
account the error bars of each correlation and the cross correlation between templates. 



Template 




Clli 




Cll2 




CI3 




CIS 




CI6 


WMAP.K 














|b 


> 30° 














Excess of 




30.0 ± 1 


D 


30.6 ± 1.7 


20 


1 ± 2.0 


IS 


4 ± 2.5 


18 


1 ± 4.9 


11.3 ± 


4 












\h 


> 40° 














Excess of 


signal 


30.2 ± 1 


9 


31.0 ± 1.8 


21 


± 2.0 


IS 


7 ± 2-8 


18 


2 ± 5.1 


11.1 ± 


4 












\b 


> 50° 














Excess of 


signal 


28.4 ± 2 





28.8 ± 1.8 


19 


5 ± 2.2 


16 


9 ± 2-9 


16 


± 4.9 


10.7 ± 


5 



4.3 Correlations with 408 MHz and 1420 MHz data 

In the upper section of Tab. |7] we list tlie correlation values be- 
tween COSMOSOMAS,WMAP and the 408 MHz template. Sig- 
nificant correlations are found at the three Galactic latitude cuts 
at all frequencies. The map at 408 MHz is regarded as a good 
tracer of synchrotron at large scales. However, at the angular scales 
left by the atmospheric filtering in the COSMOSOMAS data pro- 
cessing, the 408 MHz template includes a significant contribu- 
tion from extragalactic radiosources as can be seen from direct 
comparison with the NVSS template plotted in Fig. [16] see also 
( [Burigana et al. 2006 )l. In order to disentangle the contribution of 
extragalactic radiosources and Galactic synchrotron at high Galac- 
tic latitude, we performed a simultaneous multi-correlation analy- 
sis including the NVSS template. The results are given in the lower 
section of Tab. |7] and are interpreted as evidence for a rather wealc 
Galactic synchrotron signal left at our angular scales.. Although, 
the dispersion of the values is high for a precise determination of 
the synchrotron spectral index, the results between 11 GHz and 
22.5 GHz are compatible with a temperature spectral index of -3, 
characteristic of synchrotron emission. 

We have also studied the correlation of our maps with a des- 
ourced (Ds) Haslam et al. map given by the LAMBDA archive]^ 
The results improve in the sense that a lower correlation with NVSS 
is found, but synchrotron is still dominated by extragalactic ra- 
dio sources at our angular scales. A very similar conclusion fol- 
lows when using 1420 MHz data (Reich & Reich 1982). Finally, 
we have built a COSMO 11 Q map, combining the data from the 
two channels, and performed a cross-correlation with the Wollaben 
et al. (2006) Q map at 1420 MHz. The potential signals are at the 
level of the statistical errors in the cross-correlation and therefore 
no conclusion could be drawn for a polarized Galactic synchrotron 
signal at 1 1 GHz. An upper limit on the polarization level is set at 
2 yK (68 % C.L.). 
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Figure 16. A visual comparison between the NVSS map (top) and the 408 
MHz map (bottom) when processed according to COSMOSOMAS strategy 
and degraded to a common 1 degree angular resolution. 



this Ha map filtered according to the COSMOSOMAS data pro- 
cessing. From the table, we obtain a conversion from Rayleighs to 
K of order 40-60 ^iK/R for both channels of COSMO 1 1 . Theoreti- 
cal predictions (Dickinson et al. 2003, Bennet et al. 2003) for high 
Galactic latitude yield, assuming a mean electron temperature of 
7000-8000 K, 51-55 ^K/R at 10.9 GHz. This is compatible with 
our results. 



4.4 Correlations with Ha data 

Free-free emission is a likely contributor at our observing frequen- 
cies, especially relevant when discussing any possible dust corre- 
lated emission. At high Galactic latitudes, absorption corrections 
can be ignored at first order and we adopt the combined map of 
Ha emission provided by ( IFinkbeiner 2003b as an accurate enough 
tracer of free-free emission. Tab. [8] gives the temperature correla- 
tion coefficient between COSMOSOMAS and WMAP data and 

^ http://lambda.gsfc.nasa.gov/product/foreground/haslam_408.cfm 



4.5 Correlations with dust data 

The temperature correlations between COSMOSOMAS/WMAP 
and several dust templates are shown in Tab.|9]for various Galactic 
latitude cuts. These results are obtained with the multifit scheme 
described above. We have verified that they change by less than 
~ l^K with respect performing a simple one to one correlation 
between the COSMOSOMAS and WMAP data and the dust tem- 
plates. Within statistical errors the results are also independent on 
the templates used in the multifit. 

For |6| > 30°, the most significant correlations are obtained 
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Table 7. Temperature cross correlation between the COSMOSOMAS and WMAP maps and the 408 MHz map. The first table is for a single correlation 
analysis. The second one for the multifit case where NVSS is included in the list of templates. The reduction in the amplitude of the correlation values is due 
exclusively to the inclusion of the NVSS template in the fit. Errors are 1 a. 



Template 


Clli 


CII2 


C13 


C15 


C16 


WMAPj: 


WMAPjCa 


WMAP.Q 


WMAP.V 


WMAP.W 


6|> .^0° 


^ON .Ml 1/ 


2 3 G : U 9 


21.0 1 0.8 


119 : 9 


8.1 -; 1,2 


I <; i 2 s 


3 3 : 2 


15:02 


0.9 : 2 


, 7 : 0.2 


O.i i 0.2 


b|> 40° 


408 MHz 


25.1 ±1.0 


26.3 ± 0.9 


11.5 ±1.1 


7.6 ± 1.4 


3.6 ± 3.4 


3.9 ± 0.2 


1.7 ± 0.2 


1.1 ± 0.2 


0.9 ± 0.2 


0.3 ± 0.2 


|6|> 50° 


408 MHz 


23.1 ± 1.1 


24.0 ± 1.0 


9.5 ± 1.2 


4.3 ± 1.5 


1.9 ± 3.7 


4.4 ± 0.3 


2.4 ± 0.3 


1.6 ± 0.2 


1.3 ± 0.2 


1.0 ± 0.2 



Template 


Clli 


CII2 


C13 


€15 C16 


WMAP_K 


WMAPJCa 


WMAP.Q 


WMAP.W 


|l>|> 30° 


408 MHz 


2.8 ± 1.0 


2.8 ± 0.9 


1.1 ± 1.1 


-1.3 ± 1.4 -5.4 ± 3.3 


0.0 ± 0.3 


-0.1 ± 0.3 


-0.1 ± 0.3 


-0.4 ± 0.3 


|61> 40° 


408 MHz 


5.6 ± 1.2 


6.0 ±1.1 


1.0 ± 1.3 


-2.2 ±1.7 —6.2 ± 4.0 


0.1 ± 0.3 


-0.2 ± 0.3 


-0.2 ± 0.3 


-0.5 ± 0.3 


|I>1> 50° 


408 MHz 


4.2 ± 1.3 


4.7 ± 1.2 


-0.6 ±1.5 


-5.0 ± 1.8 -5.4 ± 4.3 


0.3 ± 0.4 


0.2 ± 0.4 


-0.1 ± 0.3 


-0.3 ± 0.4 



Table 8. Values of the temperature correlation in /iK between COSMOSOMASAVMAP and the Ha map for different cuts in Galactic latitude, 6. Errors are 1 



Template 


Clli 


CII2 


C13 


015 


C16 


WMAPJC 


WMAPJCa 


WMAP.Q 


WMAP.W 


|l>|> 30° 


Ha 


6.4 ± 0.9 


6.8 ± 0.8 


4.1 ± 0.9 


2.5 ± 1.2 


-1.7 ± 2.9 


0.9 ± 0.3 


0.3 ± 0.3 


0.0 ± 0.2 


0.0 ± 0.2 


|61> 40° 


Ha 


2.7 ± 1.0 


2.1 ± 0.9 


1.0 ± 1.2 


4.0 ± 1.5 


0.2 ± 3.5 


0.4 ± 0.3 


0.3 ± 0.3 


0.1 ± 0.3 


0.0 ± 0.3 


|l>|> 50° 


Ha 


4.0 ± 1.2 


3.2 ± 1.1 


-1.6 ±1.3 


2.7 ± 1.7 


-1.8 ± 3.9 


0.2 ± 0.3 


0.1 ± 0.3 


-0.1 ± 0.3 


-0.1 ± 0.3 



Table 9. Values of the temperature correlation between the COSMOSOMASAVMAP maps and dust maps for different Galactic latitudes. The values are 
obtained from a multifit, see text. Errors are 1 cr. 



Template 


1420 MHz 


Clli 


CII2 


C13 


C15 


C16 


WMAPJC 


WMAPJCa 


WMAP.Q 


WMAP.W 


|fc|> 30° 


AlOO 

D1RBE08 

DIRBEIO 


525.1 ± 569.1 
518.0 ± 578.1 
616.0 ± 566.1 


9.1 ± 0.9 
11.4 zn 0.9 
9.7 ± 0.9 


10.1 ± 0.8 
12.5 ± 0.8 
11.3 ± 0.8 


4.4 ± 0.9 
5.8 ± 0.9 
3.7 ± 0.9 


4.9 ±1.1 
6.3 ± 1.2 
1.9 ±1.2 


7.0 ± 2.7 
5.9 ± 2.9 
5.2 ± 2.9 


2.7 ± 0.3 

2.8 ± 0.3 
2.1 ± 0.3 


0.7 ± 0.3 
0.7 ± 0.3 
0.5 ± 0.3 


0.3 ± 0.2 
0.3 ± 0.2 
0.2 ± 0.2 


-0.1 ± 0.2 
-0.2 ± 0.2 
-0.3 ± 0.2 


|i>l> 40° 


AlOO 

DIRBE08 

DIRBEIO 


-617.0 ± 663.0 
-955.0 ± 663.0 
-314.0 ± 657.1 


6.2 ± 1.0 
6.1 ± 1.1 
4.7 ± 1.0 


7.2 ± 1.0 
7.4 ±1.0 
6.2 ± 0.9 


0.4 ±1.2 
1.0 ±1.2 
1.4 ±1.2 


3.4 ± 1.4 
2.3 ±1.4 
-0.7 ± 1.5 


1.5 ± 3.5 
0.0 ± 3.5 
-0.6 ± 3.5 


1.5 ± 0.3 
1.2 ± 0.3 
0.8 ± 0.3 


0.6 ± 0.3 

0.5 ± 0.3 
0.3 ± 0.3 


0.3 ± 0.3 
0.2 ± 0.3 
0.2 ± 0.3 


0.1 ± 0.3 
0.0 ± 0.3 
-0.1 ± 0.3 


|6|> 50° 


AlOO 

D1RBE08 

DIRBEIO 


-1487.0 ± 732.1 
-1660.0 ± 731.1 

-651.1 ± 723.0 


2.6 ± 1.2 
1.6 ± 1.2 
2.8 ±1.1 


1.8 ±1.1 
1.6 ± 1.1 
4.4 ±1.0 


2.6 ± 1.3 
2.6 ± 1.3 
3.6 ± 1.3 


3.6 ± 1.6 
2.0 ±1.6 
0.2 ±1.6 


-2.8 ± 4.0 
-5.0 ± 4.0 
-5.0 ± 3.9 


1.4 ± 0.3 
0.9 ± 0.3 
0.8 ± 0.3 


0.5 ± 0.3 
0.4 ± 0.3 
0.2 ± 0.3 


0.3 ± 0.3 
0.2 ± 0.3 
0.2 ± 0.3 


0.2 ± 0.3 
0.0 ± 0.3 
-0.1 ± 0.3 
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for both COSMO 11 channels, with a signal to noise ratio greater 
than 10, and for the WMAP K channel, with a signal to noise 
around 7-9. Marginal detections are also found with the Ka chan- 
nel. No significant detection is found at the higher frequency chan- 
nels of WMAP. The channels of COSMO 15 also exhibit some clear 
detections, although for some templates these are less significant. 
At higher Galactic latitudes the amplitude of the coiTelated signal 
decreases and for j&j > 50° the detections are below the 3(j level. 
Overall, we find clear evidence for a dust correlated high Galactic 
latitude microwave emission in the frequency range 1 1-33 GHz. 

We note that the two COSMO 1 1 channels present a small sys- 
tematic difference in the values of the correlations at > 30° and 
|6| > 40°, of order 1-1.5 /iK, similar to the statistical error. In prin- 
ciple, this could be due to a significant level of polarization in the 
dust-correlated signal but more sensitive maps are required in order 
to establish whether this is a reliable difference. 

The dust-correlated signal detected in Tab. |9] presents a spec- 
tral index between 11 and 33 GHz of -2.1 ±0.2, incompatible with 
synchrotron emission (~ —3) but consistent with free-free emis- 
sion. However, the remarkable flattening found in the frequency 
range of the COSMOSOMAS experiment cannot be understood in 
terms of free-free emission properly, suggesting an improper trac- 
ing of the free-free emission at high Galactic latitude by the Ha 
template used in the multifit analysis. This would be the case in re- 
gions of high dust emission where dust also causes high extinction 
in the visual. A careful study of the regions that mostly contribute 
to the dust-correlated signal in the COSMO 1 1 data at |6| > 30° re- 
veals areas with a very strong emission in the DIRBE maps which 
indeed are not Hq emitters. Dust extinction in these regions would 
cause that free-free emission would not be traced by the Ha tem- 
plate affecting to the correlation values listed in Tab. [9] In order to 
explore if this is the case, we identified regions with very strong 
dust emission at high Galactic latitudes and performed a new cor- 
relation study by extending the mask to include those pixels whith 
a significant correlated emission with the COSMOSOMAS maps. 
For definiteness, we extended the mask to those pixels where the 
amplitude of the normalised covariance between the channel 2 of 
COSMOSOMAS (the one with the best sensitivity) and DIRBE08 
(i.e. < tc2tD08 > /(cr(tc2)cr(tD08))) was greater than 5. Adjacent 
pixels were also masked. In total, only a few regions of size 1 — 5 
square degrees were removed from the analysis and a new multifit 
correlation analysis was performed The results are listed in Tab.|10| 
We find a decrease in the amplitude of the dust- correlated signal at 
any frequency with respect to the case in Tab.|9]but the correlations 
at > 30° are very significant for COSMO 1 1 and WMAP_K and 
significant for the other COSMOSOMAS channels. There is also 
a trend to lower amplitudes as we move to high Galactic latitude. 
In order to see if this dust-correlated microwave emission could be 
associated with synchrotron, free-free or a different emission pro- 
cess we conducted a detailed analysis of the DIRBE08 correlation 
161 
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Jfor 16| > 30°. 

A maximum likelihood approach has been used to find the 
best-fit parameters within three different models. First, we con- 
sidered a single power law with two free parameters (,f{u) — 
Aiy /22 GHz)*') as a fitting function; second, we considered a two- 
component fit, in which one component is given by a power-law, 
and the other is described by the standard combination of Draine 
& Lazarian models (CNM, WNM and WIM, with relative ampli- 
tudes of 0.43, 0.43 and 0.14, respectively). Finally, we considered 
a phenomenological approach and tried to fit the data with a sin- 
gle Gaussian law given by f{v) — A exp {—{v — v^f' j (2it^)) in 
order to mimic the behaviour of a general spinning-dust model. 
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Figure 17. Correlation between COSMOSOMAS. WMAP channels and 
the DIRBE08 map in flux density units (Jy/sr) for |6| > 30° after masking 
some localized regions, see text for details. Oveiplotted is the best fit model 
of a Gaussian model with parameters A = 23.6Jy/sr, t'o = 21.7 GHz and 
o- = 15.8 GHz. 



The best fit is obtained in the third case, with a reduced value 
of x^/dof = 0.72 (probability of 63%), although the other two 
cases show similar values for the goodness-of-fit, obtaining 1.05 
(39.1%) and 0.93 (46.9%) for model 1 and 2, respectively. We note 
that in the case of model 1 the index of the power-law (p = 0.1 ± 
0.2)) rules out synchrotron emission but is compatible with free- 
free. 

For the third model, the best fit parameters (obtained from the 
marginalized likelihood function) are A = 23.5^3 ^ Jy/sr, fo ~ 
21.713.7 GHz and a = 15.813;° GHz. In Fig.[l7]we plot the data 
of Tab.fTOland this model. 

This "blind" best-fit function peaks at a slightly lower fre- 
quency than the considered Draine & Lazarian model, and has a 
slightly broader width. In the case of interpreting this anomalous 
signal in terms of a spinning dust model, this would have an im- 
pact on the actual size distribution of the grains responsible for the 
emission. 

In summary, these results appear to support the existence of 
a spinning dust microwave emission process at high Galactic lati- 
tudes. In order to disentangle whether this anomalous microwave 
emission is associated to well defined regions as it is the case of 
the Perseus molecular complex (Watson et al. 2005), or to a diffuse 
large-scale Galactic component, maps of higher sensitivity and res- 
olution are required in the frequency range 10-30 GHz. 



5 CONCLUSIONS 

We have presented new results of the COSMOSOMAS experi- 
ment at 1 1 GHz for a sky coverage of ca. 6500 sq. deg. A cross- 
correlation analysis using this new dataset, previous COSMOSO- 
MAS data in the frequency range 12-17 GHz, the 3rd year WMAP 
data and several foreground templates yields the following results 
at high Galactic latitude (jbj > 30°): 

• The presence of a CMB component of amplitude 27±2 /^K 
in the COSMOSOMAS channels in agreement with the expected 
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Table 10. Temperature correlations between the COSMOSOMAS and WMAP maps and dust maps for different Galactic latitudes. The values are obtained 
extending the mask to the particular regions where there is a high correlation between COSMOSOMAS and DIRBEOS, see text. Errors are 1 a. 



Template 


1420 MHz 


CI 


02 


CI3 


CIS 


CI7 


WMAPJt 


WMAP.Ka 


WMAP.Q 


WMAP.W 




b|> 30° 


AlOO 

DIRBEOS 

DIRBEIO 


96.1 ± 571.0 
-139.1 ± 573.1 
115.0 ± 568.1 


4.7 ± 0.9 
4.6 ± 0.9 

4.8 ± 0.9 


5.9 ± 0.8 
5.8 ± 0.8 
6.6 ± 0.8 


2.4 ± 0.9 
2.8 ± 0.9 
1.2 ± 0.9 


2.3 ±1.1 
3.5 ±1.2 
-0.6 ± 1.2 


4.5 ± 2.7 
2.5 ± 2.9 
2.8 ± 2.9 


1.9 ± 0.3 
1.8 ± 0.3 
1.3 ± 0.3 


0.4 ± 0.3 
0.4 ± 0.3 
0.3 ± 0.3 


0.2 ± 0.2 
0.2 ± 0.2 
0.1 ± 0.2 


-0.1 ± 0.2 
-0.2 ± 0.2 
-0.3 ± 0.2 


|6|> 40° 


AlOO 

DIRBEOS 

DIRBEIO 


-962.0 ± 663.1 
-1252.0 ± 663.0 
-443.1 ± 659.0 


3.5 ± 1.0 

3.4 ±1.1 

3.5 ± 1.0 


3.7 ± 0.9 
3.9 ± 1.0 
4.7 ± 0.9 


-0.4 ± 1.2 
0.3 ± 1.2 
1.1 ± 1.2 


1.9 ± 1.5 
0.8 ± 1.5 
-1.4 ± 1.5 


0.6 ± 3.5 
-1.2 ± 3.5 
-1.0 ± 3.6 


1.4 ± 0.3 
1.0 ± 0.3 
0.7 ± 0.3 


0.6 ± 0.3 
0.4 ± 0.3 
0.2 ± 0.3 


0.3 ± 0.3 
0.2 ± 0.3 
0.1 ± 0.3 


0.2 ± 0.3 
0.1 ± 0.3 
-0.1 ± 0.3 


|l>|> 50° 


AlOO 

DIRBEOS 

DIRBEIO 


-1541.0 ± 733.0 
-1710.0 ± 732.0 

-656.1 ± 723.1 


1.8 ± 1.2 
0.7 ± 1.2 
2.7 ±1.1 


0.7 ±1.1 
0.4 ±1.1 
4.3 ± 1.0 


2.5 ± 1.3 

2.5 ± 1.3 

3.6 ± 1.3 


3.5 ± 1.6 

2.0 ± 1.6 
0.2 ± 1.6 


-2.8 ± 4.0 
-4.9 ± 4.0 
-4.8 ± 3.9 


1.3 ± 0.3 
0.9 ± 0.3 
0.8 ± 0.3 


0.5 ± 0.3 
0.4 ± 0.3 
0.1 ± 0.3 


0.3 ± 0.3 
0.2 ± 0.3 
0.2 ± 0.3 


0.2 ± 0.3 
0.0 ± 0.3 
-0.1 ± 0.3 



CMB fluctuation after removal of the Fourier terms implied in the 
atmospheric filtering. 

• For the angular scales of our experiment, unresolved extra- 
galactic sources are found to be the dominant foreground at 11 
GHz. The measured contributions of unresolved radiosources in the 
range 1 1-23 GHz are in very good agreement with predictions by 
de Zotti et al. (2005). 

• The synchrotron component at high Galactic latitude is found 
with an amplitude of 3-6 /xK at 1 1 GHz, and decreases at higher 
frequency with a temperature spectral index compatible with -3. 

• The cross-correlation with the H^ map gives an amplitude of 
3-6 /iK at 11 GHz. Correlations with higher frequency COSMO- 
SOMAS and WMAP data verify a temperature spectral index con- 
sistent with -2, characteristic of free-free emission. 

• A dust correlated emission is detected in each of the COSMO- 
SOMAS channels at |6| > 30°. The amplitude of the signal ranges 
from 10-12 /xK at 11 GHz down to 4-7 fj,K in the 12-17 GHz and 
2.1-2.8 /iK at 22 GHz. The Galactic latitude dependence supports 
a Galactic origin for this signal. An important fraction of this cor- 
related signal at 1 1 GHz comes from regions of high dust emission 
where free-free emission is not well traced by the Hq template due 
to extinction. After masking those regions the remaining dust corre- 
lated signal detectable in the frequency range 1 1-33 GHz - of order 
5-6 /iK at 1 1 GHz - shows a clear flattening, which is not com- 
patible with the classical spectral index of synchrotron emission. 
This correlated signal can be described by models resembling spin- 
ning dust emission, with an associated flux density peaking around 
20 GHz. 



ACKNOWLEDGEMENTS 

We are thankful to L. Toffolatti and J. Gonzalcz-Bucno for very 
useful discussions on unresolved radiosource contributions at the 
frequencies of COSMOSOMAS. We also thank J.L. Salazar and 
the personnel of the LAC engineering division, instrument main- 
tenance and operations of Teide Observatory for their support to 
the COSMOSOMAS experiment. This research has been partially 
funded by project AYA 2001-1657 and AYA 2005-06453 of the 
Spanish Ministry of Education and Science. 



16 



S. R. Hildebrandt et al. 




Figure 18. COSMOll scanning strategy together with the two polarization directions observed in Galactic coordinates. The image used for the sky is that 
provided by WMAP for its lowest frequency channel (K-band, 22.5 GHz), which is the closest to COSMOll observed frequencies. Top: image centered at 
Galactic longitude 0° . The Celestial north pole is again marked with an "N". Cygnus A complex is labelled with a "C". HEALPix pixelisation scheme is used, 
IGorski et al. 2005> . 



COSMOSOMAS 10-12 GHz Data. 
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